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Input: pts : vector of points, n : number of points, Of ,d 9 O : parameters and std-dev of fitted line. 
Output : true if curve contains a zig-zag- shape, false otherwise. 
Sequence : 

if a > thresho then 

valO = pts [0] .x* cos or +pts [0] .y* sin a + d 
count = 1 
for i = 1 to (n-1) do 

val =pts[i].x* cos or +pts[i].y* sin Qf+c/ 
if(ra/*ra/0 <=0)then 
ra/0 = va7 

COZJW^ = 1 

else 

count = count + 1 

if {count >= min-points-for-zig-zag-shape ) then 

return true 
endif 
endif 
endfor 
endif 

return false 
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Algorithm FindSeedRegion 

Input: lines [i]: vector of lines for esch image row (or column) /, 

n : number of image rows (or columns) 
Output : set of lines (seed region) or empty set (no seed found). 
Sequence : 
for i = 0 to (n-3) do 

for lj in lines fij do 

for 1 2 in lines[i+l] do 

for 1 3 in lines[i+2] do 

if overlap( l l9 l 2 ) and overlap( 1 2 ,1 3 ) then 
(n, d) = GtPlan e(l t ,1 2 ,1 3 ) 
if rms(l } ,1 2 ,1 3 )< thresh 1^ then 
seed= {l l3 l 2 J 3 } 
removeQj ,1 2 ,1 3 ) 
return seed 
endif 

endif 
endfor 
endfor 
endfor 

endfor 
return {} 
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Algorithm RegionGrowing 

Input: region : set of lines as seed region, 

lines fij: vector of lines for esch image row (or column) i, 

n : number of image rows (or columns) 
Sequence : 

A=0,b = 0 

for 1 in region do (A,b) = add(A,b,l) endfor 
(n,d) =solve(A,b) 
open = region 
while not empty(open) do 

1 7 = select(open),open = open - {1 } } 
for i in neighbor(index(l I )) do 
for 1 2 in JinesfiJ do 

if overlap(Jj ,1 2 ) and rms(l 2 ) < tbresh2 ims then 
region = region + {1 2 }, 
(A,b)=add(A,b y l 2 ), fad) = solve(A,b) 
open = open + {1 2 }, 
remove(l 2 ) 
endif 
endfor 
endfor 

endfor 

plane = {n,d,A,b,region} 
planes = planes + {plane} 
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